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Abstract 

While agreeing with our exact expression for the time dependence of the motion of a free particle 
in an initial superposition state, corresponding to two identical Gaussians separated by a distance 
d, at temperature T, Gobert et al., in the preceding Comment [Phys. Rev. A xxx], dispute our 
conclusions on decoherence time scales. However, the parameters they used to generate their 
figures are outside the regime of validity of our interpretation of the results and, moreover, are not 
of physical interest in that they correspond to T ~ 0. The point is that in their figures they have 
chosen the thermal de Broglie wavelength = hy/mkT to be equal to slit spacing d, whereas we 
have clearly stated [in the paragraph preceding Eq. (21) of our paper] that decoherence occurs and 
that our expression for the decoherence time applies only in the limit where d is large compared 
not only with the slit width a but also with the thermal de Broglie wavelength, d 3> Xth, 
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For a free particle in an initial superposition state, corresponding to a pair of Gaussians 
each of width a and separated by a distance d, at temperature T, we presented an exact 
expression for the coordinate probability, P(x, t) . Gobert et al. (hereafter referred to 
as GDA), in the preceding Comment |3[, agree with this result, which is valid for arbitrary 
temperature T and arbitrary strength of the coupling to the heat bath as measured by the 
dissipative decay rate 7 . GDA also agree that their results for the corresponding reduced 
density matrix are consistent with_our results for W(q,p,t), the Wigner function for the 



probability in Wigner phase space jj| (we note that the Wigner distribution function is the 
Fourier transform of the density matrix and contains the same information 5}). Based on 
our result for P(x, t) we concluded that decoherence without dissipation can occur for all 
cases of interest. GDA dispute this conclusion. Our purpose here is to show their analysis is 
flawed due to a choice of parameters that are not of physical interest in that they correspond 
to T pa 0. The point is that their figures correspond to a slit spacing d equal to the thermal 
de Broglie wavelength, 

But we have clearly stated [in the paragraph preceding Eq. (21) of our paper] that deco- 
herence occurs and that our expression for the decoherence time applies only in the limit 
where d is large compared not only with the slit width a but also the thermal de Broglie 
wavelength, d 3> Xth, Indeed, decoherence is fully developed even for the very modest 
ratio d/Ath- To illustrate this we show in Fig. 1 a plot of P(x,t) [multiplied by o to make 
it dimensionless] versus x (divided by a) at a time t = 2mad/5h = t m ix/5 for three different 
values of the ratio d/X t ^. The solid curve corresponds to d/A t h = 5 and there we see that 
there is no hint of an interference pattern, that is, decoherence has occured. The two dashed 
curves, which are nearly indistinguishable, correspond to the GDA choice d/A t h = 1 and 
the zero-temperature case d/X t h = 0, the zero temperature case having the slightly larger 
interference amplitude. In all three cases we have neglected the coupling to the bath (that 
is, set 7 = 0) and made the GDA choice d = 20a. We assert that this Figure refutes those 
of GDA, showing that under the conditions we have stated decoherence without dissipation 
does occur. GDA ask how can it be that the interference pattern has not disappeared in their 
Fig.l. The answer of course is that they have chosen values of the parameters corresponding 
to an effectively zero temperature, for which in the absence of dissipation there is of course 
no decoherence. As our Fig. 1 shows, if they were to choose only a modestly higher tern- 



perature, corresponding to dj A t h = 5, they would see no interference. We forbear to present 
a repeat of GDA's Fig. 1, since for zero temperature the figure would be indistinguishable 
from theirs, while for <i/A t h = 5 the figure would be utterly featureless, consisting of two 
Gaussians propagating independently without interference. Perhaps a numerical illustration 
might be helpful. Since the magnitude of the de Broglie wavelength A^ plays a key role in 
the analysis, we use (JTJ) to write 



As an illustrative example, Zurek, in his oft-cited paper jfj, chooses m = lgm, d = lcm and 
T = 300K, corresponding to the very large ratio d/A t h = 2 x 10 20 . Put another way, the 
GDA choice of d/Xth = d would require a temperature of 8.1 x 10 _39 iT! This is perhaps an 
extreme choice of parameters, but it illustrates the fact that the condition d ^> A^ is not 
unique to our discussion but has generally been understood to apply in earlier discussions 
of decoherence. In general, the decoherence times are typically much smaller than the 
decay time 7 _1 . But this is so only under our condition Ath d, no matter how one defines 
decoherence. Next, we turn to the question of how decoherence should be defined. As we 
have stressed more than once, "-a quantitative measure of decoherence depends not only 
on the specific system being studied but also on whether one is considering coordinate, 
momentum, or phase space" |4|. On the other hand, based on their density matrix analysis 
(which is nothing more than that given previously by us for Wigner phase space [4] ), GDA 
claim that this implies that there is no decoherence in physical space, which as we have 
seen is not true. Moreover, while again recognizing that there are various ways of defining 
decoherence, our definition of an attenuation coefficient a(t) was such as to ensure that 
a(t) = 1 for all times in the case where T = and 7 = 0, corresponding; to a pure Schr 
odinger cat state at all times (as discussed in detail in Sec. VI B. of |7[). This choice 
leads to a time scale which is convenient for this system but other definitions of decoherence 
time (which will also involve the temperature) might also be considered, depending on what 
particular features of the time development are of interest. For example, it is instructive to 
write our general result for P(x,t) in the form Q, 0] 






(3) 
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where Pq is the probability distribution for a single wavepacket, given by 

p ° (x ' t)= ^W) exp {-A)}- (4) 

Here and in Q w 2 (t) is the variance of a single wavepacket, which in general is given by 

w\t)=a 2 - [m ^ t)]2 +s(t), (5) 
where a 2 is the initial variance, [x(0),x(t)] is the commutator, and 

s(t) = ({x(t) - x(0)} 2 ) , (6) 



is the mean square displacement. Also, a(t), which can be defined as the ratio of the factor 
multiplying the cosine in the interference term to twice the geometric mean of the first two 
terms 2j is given by the following exact general formula 

n 

As we emphasized in |2(, it is only for the case of Ohmic dissipation, high temperature 
(kT >> Hj) and d » Xth, that (J2J) reduces to 

t- 2> 



a(t)^exp|-^-j j, (8) 
where 

Td = 1 — , (9) 

'kT/m d 



the result quoted by GDA in their equation (4). But, this result clearly does not apply for 
d = Xth, the GDA choice of parameters. 

GDA continually refer to FLO's "- imperfect preparation of the initial state -". What 
perhaps they mean is that the initial state of FLO is not a pure state, which is certainly the 
case. But a pure state necessarily corresponds to a particle at zero temperature. FLO argue 
that the picture of a particle at zero temperature suddenly coupled to a heat bath at a high 
temperature T is not a realistic picture of the physical situation. The essence of FLO's paper 
is to show that one can, by measurement, prepare an initial state that is entangled with 
the bath, with a temperature equal to that of the bath, but with a probability distribution 
identical with that of the pure state. They argue that this is a more realistic initial state. 
We conclude that the comment of GDA is simply misleading due to their unrealistic choice 
of parameters and we reiterate our previous conclusion that decoherence can occur without 
dissipation. 
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FIG. 1: P(x,t) the coordinate probability at time t (multiplied by a to make it dimensionless) as a 
function of x (divided by a) is shown at a time t = t m \ x /h. The solid curve, which shows no hint of 
interference, corresponds to At = d/5. The two dashed curves, which are nearly indistinguishable, 
correspond to the GDA choice of parameters: At = d and to the zero temperature case At = oo, 
with the zero temperature interference pattern having a slightly larger amplitude. 
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